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Abstract. - We have measured the mass diffusion coefficient D of ^He in superfluid *He at 
temperatures lower than were previously possible. The experimental technique utilizes scin- 
tillation light produced when neutrons react with "^He nuclei, and allows measurement of the 
^He density integrated along the trajectory of a well-defined neutron beam. By measuring the 
change in '^He density near a heater as a function of applied heat current, we are able to infer 
values for D with 20% accuracy. At temperatures below 0.7 K and for concentrations of order 
10~*, we find D = [2.0^^ 2] X'-(6-5=Fi-2) gjn^/s, in agreement with a theoretical approximation. 



Lntroduction. - We have developed a new technique for studying the distribution of 
■^He atoms at low concentrations in superfluid ^He and have appUed it to the measurement 
of the '^He mass diffusion coefficient D at temperatures below 1 K. Previous measurements 
of D have not extended to temperatures below 1 K The earlier results are based on 

thermal conductivity measurements and, because D becomes very large at low temperatures, 
an impractical amount of heat input is required to achieve a measurable temperature difference 
below 1 K. 

We are interested in the magnitude of D near and below 0.5 K, where it is expected to 
be large Reference 0] indicates that D — .025 cm^/s at T = 1.2 K, and that it increases 
rapidly with decreasing temperature; the temperature dependence of D is predicted to be T"^ 
for T < 0.6 K ||. 

Theory of the experimental technique. - Our method of measuring D is based on a 
comparison of the "^He distribution in a cell with and without a source of heat in the liquid. In 
the treatment below we consider heat flow and ^He concentrations that are sufficiently small 
such that we can assume the liquid is isothermal. We also ignore the relatively infrequent 
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^He-'^He collisions that do not influence the steady-state "^He atom distribution imposed by 
the heat flow. This system can be analyzed in terms of the usual macroscopic two-fluid model 
for liquid ^He below the A point Q ; these two fluids comprise the superfluid component that 
carries no entropy, and the normal component that, at sufficiently low temperatures, carries 
thermal energy as a pure phonon gas. When a heat source and a heat sink are introduced, 
the heat flux q associated with the normal component flow is 

q = psTvn (1) 

where p is the "^He liquid density, s is the entropy per unit mass, T is the temperature, and 
Vn is the normal component velocity. The "^Hc atoms scatter from the phonons and are thus 
carried along with the normal fluid flow. Consequently, the '^He atoms accumulate at the heat 
sink. The '^He atoms also scatter from the background phonons that exist when the superfluid 
'^He is at finite temperature; this scattering tends to randomize the "^He distribution. The net 
effect is described by the diffusion equation 

Xv„ - DVX = (2) 

where X represents the spatially varying ^He concentration (see, e.g., §29, eq. (13) of [Q). In 
steady state, 

^ - = -V • /7.sv„. (3) 

We assume p and s are spatially uniform under isothermal conditions in which case V • Vn = 0. 
Furthermore, if the flow is laminar, then V x Vn = 0, which allows us to write Vn = V$, 
where the velocity potential $ satisfies Laplace's equation, = 0. Substituting this result 
into eq. (|), 

XV$ - DVX = ^ = ^ VX = V log(X), (4) 

which has the solution 

X = Ce*/^. (5) 

The coefficient C, determined by the requirement that the total number of '^He atoms must 
be constant for a closed system, depends on the total heat input and D. For the case of a 
point heat source of power P in a large bath (i.e., large enough so we can ignore boundary 
conditions), 

( ) ^ =^ ( ) '^^^'^ ^ ) P ^ /Q\ 

Anr^ psT psT Aiir Anr 

where r is the distance from the source, and a = P/{psT). In this case. 



X(r) = Ce-"'/^""'' (7) 

and D can be determined by measuring the spatial extent of the '^He-depleted region as a 
function of P. Alternatively, D can be inferred if the integrated '^He density along a specific 
path near the heat source is measured as a function of P. The latter approach serves as the 
basis of the measurements reported here. 

Description of the experiment. ~ Use of a neutron beam to determine tomographically the 
distribution of ^He atoms subjected to heat currents in a superfluid-^He-filled cell is described 
in Appendix C of For the measurements reported here, we use a variant of this technique; 
we determine the differential heat on/off effect using a fixed and well-characterized neutron 
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beam. A scale drawing of our apparatus is shown in Fig. ^ The 100 cm'^ ceh containing 
superfluid ^He is mounted on the end of a horizontal dilution refrigerator; the OFHC copper 
mixing chamber of the dilution refrigerator, which serves as a heat sink, forms one end of the 
cell. The cell body is a cylindrical shell of aluminum approximately 1 mm thick. A fused 
silica window forms the other end of the cell. 

Neutrons are absorbed by '^He according to the reaction 

^He + n ^ p + t + 784keV (8) 

and the energy released to the proton and triton create XUV (80 nm) scintillation light as 
they are stopped in the superfluid '*He. The XUV light is wavelength-shifted to 450 nm with 
tetraphenyl butadiene (TPB) evaporated onto the cell walls and onto an aluminized mylar 
sheet that covers most of the copper end of the cell. The 450 nm light is transmitted through a 
series of windows. In addition to the cell window, there are two fused silica windows at about 1 
K, one fused silica window at about 77 K, an acrylic (blackbody radiation absorbing) window 
at 77 K, and a fused silica window at room temperature. The scintillations are then detected 
by a Hamamatsu R329-02 photomultiplier, set to a threshold of 5 ± 0.5 photoelectrons; the 
detection efficiency of "^He capture events varies spatially within the cell in the 40 to 100% 
range, with the variation depending primarily on the distance from the mylar sheet. 

The experiment was performed at the Los Alamos Neutron Science Center spallation source 
and used Flight Path 11 A, which comprises a cold neutron guide that views a liquid hydrogen 
moderator. The neutron beam was formed by coUimating the output of the cold neutron guide 
to a 2.5 mm FWHM beam with a divergence less than 1 x 10""^ rad. The pulse structure of 
the source is not crucial to the technique, but was used in our measurements to eliminate the 
background due to the intense gamma and energetic charged particle pulse associated with 
the proton current pulse incident on the spallation target. The pulse counters were gated open 
in a time window of 6 ms following a delay of 17 ms from the spallation pulse; the distance 
between the experiment and moderator was 24 m, so the neutron velocity range was 1040 to 
1410 m/s. Approximately 3% of the incident neutrons are scattered by the superfluid '*He, 
and at our level of accuracy, the effects due to scattered incident neutrons can be neglected. 

The initial '^He concentration Xq was varied during the course of the experiment such that 
between 5 and 20% of the incident neutrons were absorbed, corresponding to 7 x 10~^ < 
Xq < 3 X 10^"*. The gated count rate was typically 450/s {Xq = 10^"*) with a background 
of 12/s. The background was determined by operating the experiment with the cold neutron 
component of the beam blocked by thin Cd and ^Li doped plastic sheets, and was independent 
of heat input and relatively insensitive to the '^He concentration. This background was largely 
due to long-lived molecular optical excitations in the superfluid bath j|] that were left over 
from the spallation radiation flash. In addition, the neutron beam intensity transmitted 
through the cell was measured with a ^Li-doped glass scintillator detector, gated identically 
with the helium scintillation signal. The proton pulse current provided by the facility was 
used to obtain a relative normalization of the cold neutron intensity. 

In order to perform tomography, the cryostat containing the dilution refrigerator and cell 
was placed on a motor-controlled table that allowed horizontal and vertical motion perpen- 
dicular to the neutron beam. Linear resistive displacement transducers linked to the table 
motion allowed positioning to 0.5 mm absolute accuracy, and to better than 0.25 mm repro- 
ducibility. The heat source was a carbon composition resistor (50 at room temperature) 
located near the bottom of the cell and supported by an insulating alumina column. A cali- 
brated ruthenium oxide temperature-sensing resistor supported by its wire leads was located 
near the mixing chamber end of the cell as shown in Fig. ^ This resistor was used in con- 
junction with a heater located in the mixing chamber to provide active relative temperature 
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stabilization to within 20 niK under varying heat input conditions. The uncertainty in the 
absolute temperature was about 0.03 K and depended on the operating temperature. The 
maximum power delivered to the cell was in the 7 to 15 mW range and was limited by the 
temperature-dependent cooling power of the refrigerator. 

The helium scintillation signal provides a relative measure of the volume-integrated "^He 
density along the neutron beam. In order to interpret the effects of the applied heat flux it 
is necessary to determine the velocity potential ^{x,y, z) within the cell. This was done by 
numerically solving = on a 100 x 100 x 100 grid using standard relaxation techniques, 
with the boundary condition that the perpendicular component of Vn was zero at all bound- 
aries, except at the mixing chamber, where the condition $ — was set. The perpendicular 
component of Vn at the resistor surface was determined by the total heat input and assumed 
to be uniform. Finally, the normalization was determined by integrating e* over the grid, 
with D left as a parameter (scale factor) to be determined experimentally. The integrated 
^He density was determined by evaluating the integral of X over the neutron, as beam shown 
in Fig. |l|. 

We verified experimentally that eq. does indeed describe the distribution of '^He in 
the cell subject to heat flow (to be published). For determining D, the beam was directed 
through the cell in the region near the heater; the beam location used for the measurements is 
indicated in Fig. |^. Results of the numerical calculation of the integrated ^He density along 
the neutron path shown in Fig. Q are given in Fig. ^ as a function of heat power; the value of 
JD at a given temperature enters as an unknown scale factor on the horizontal axis of Fig. ^. 
By measuring the relative rate of change of scintillation rate R (which is proportional to the 
relative change in integrated ^He density) near zero power, D can be determined from 

-0.078 cm-i 

" dR/da ' ^> 

where the numerical factor is the maximum slope (i.e., at zero heat power) in Fig. ^. Because 
we expect that D scales as at low temperatures (see below), we can cast the data in terms 
of a universal function if we plot the experimentally-observed relative change in scintillation 
rate as a function of x = aT^. Then eq. (^, with a replaced by x, yields Dt directly, where 
D = DtT-\ 

Experimental Results. - The data set is shown in Fig. ^. The relative scintillation rate 
R{x)/R{0) as a function of heater power (parameterized by x) was determined by subtracting 
the background rate from the observed scintillation rate, and dividing by the background- 
subtracted scintillation rate at zero power. The vertical errors were determined by counting 
statistics, while the horizontal errors were determined by the estimated uncertainties in T (and 
therefore s). The low temperature points tend to fall on the same curve; the small values of 
at the lower temperatures, along with the limited cooling power of the dilution refrigerator, 
restricted these data to a; < 2.2 as shown in Fig. ||. This is sufficiently small so that the 
variation in relative scintillation rate can be assumed to be linear over the range < a; < 2.2 
(see Fig. |^). The low temperature points (e.g., the data below 0.7 K) are re-plotted in Fig. ^ 
using expanded axes. The points fall on a single line with reduced = 0.78. From the slope 
of this line, along with eq. (||), we arrive at 

D = DtT-'' = {1.6 ±0.2)T-''. (10) 



As expected, there was no observed effect due to concentration (e.g., the low temperature 
data are consistent as indicated by x^), so we present no analysis on the effect of varying Xq. 
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Also, from eq. (|l|), u„ < 10 cm/s for all the data; this is sufficiently small so that turbulence 
effects are unlikely to be important. 

Finally, D was determined at each temperature independently as shown in Fig. ^, along 
with results from previous work. It can be seen that our results extrapolate to the Beenakker 
et al. [0 results, but both experiments are inconsistent with the Ptukha ||] results. One 
can further analyze the lowest-temperature data in a two-dimensional fit using D = DtT~''^ 
yielding 

[2.0+^-^1 T-(6.5T1.2) ^^^2/ .^1) 
[ -I.2J ' ^ ' 

and indicates the level at which the theory presented below is constrained by the experiment. 

Microscopic theory of ^ He impurity diffusion. - The expectation that D will vary as 
T^^ when collisions with phonons provide the dominant '^He scattering mechanism (T < 0.6 
K) was briefly discussed in Ref. [|| (Sec. 7.2); here we add further details to that argument. 
From Fig. 2 of we infer that the '^He-phonon collision frequency l/rph varies as T"'. The 
data at T = 0.5 K and concentration X = 1 x in this figure imply 

1/Tp,, = 4.8 X 10^" X T'' s~^ (12) 

Because the number of collisions per unit volume per unit time must be the same for phonons 
and ^He atoms, we have 

Tph T3 "3 Tph 

where rt3 = 2 x 10^^ X cm^'^. The phonon density is approximately (e.g., Q, §16, eq. (28)) 

uph = 8nC{3){kT/hcf = 2 x lO^^T^ cm^^ (14) 

where C(3) = 1.202, c — 2.4 x 10^ cm/s is the phonon phase velocity, k is Boltzmann's constant, 
and h is Planck's constant; therefore 

— = 4.8 x lO^r^s"^ (15) 
T3 

The collisions measured in Ref. are those effective in phonon transport, i.e., those collisions 
that change the phonon momentum by a significant amount. For elastic collisions, which 
dominate here, the momentum transfer for a phonon of momentum q scattered through an 
angle 6q (where Aq = 2qsm9q/2), is equal to the momentum transfer to a '^He atom with 
momentum p3, scattered through an angle 63 (where Ap3 = 2p3sin03/2), in a given collision. 
For each collision, 

pI{1 - cos 03) - q\l - cos0q) ^{1- cos 03) « (1 - cosdg)(q^) / {pi); (16) 

where (...) indicate a thermal average. From elementary considerations, the average squared 
momenta are 

„.),(^,iM^I)!, (pi) = «T™; (17) 

where = Iflm^ is the effective ■^He mass in superfluid ^He (7713 is the ^He atomic mass). 
Numerically, {q^) / {pf) = 7.5 x lO'^T so for T < 1 K, pi » q^ and (cos 61,) « 0. From 
kinetic theory, the diffusion coefficient is proportional to the product of 1/(1 — cos^s) and 
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(^'3^3) = T3{vl), where V3 is the ^He velocity and A3 — is the distance a ^He atom travels 
between collisions: 

D = J(A3«3)/(1 - cosOg) = i(t|)T3/(l - cos^,). (18) 

We find that 

D ^1.2 T-'^cm^/s^ DtT^'^. (19) 

The factor Dt is likely correct to within an order of magnitude and is applicable to temper- 
atures near and below 0.6 K, e.g., where the thermodynamic properties of superfluid "^He are 
dominated by phonons. More importantly, however, this result indicates the functional form 
of D at low temperature. The value of Dt from eq. ( |l9|) is in reasonable agreement with 
our experimental result, eq. (|lo|). A useful figure-of- merit for the effectiveness of a given heat 
input power for transporting '^He atoms by phonon-dominated entrainmcnt is, by use of eqs. 
(|) and (0), sT/D oc (s cx T'^ for T < 0.7 K). The deviations from the low temperature 
T^^ behavior of D that are evident for T > 0.65 K in Fig. || are likely due to the influence of 
of rotons at higher temperatures. 

Conclusion. - We have employed a novel differential absorption technique to determine 
the mass diffusion coefficient of '^He in superfluid ^He as a function of temperature. The 
measurements extend to much lower temperature than in any previous work. We have verified 
that the diffusion coefficient is closely proportional to at temperatures below 0.7 K. 

This work was supported by Los Alamos National Laboratory LDRD under Project Num- 
bers 97041, 98039, and 2001526DR. 
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Fig. 1 - Scale drawing of the test cell. There are indium seals at the connections between the aluminum 
cell cylinder, demountable window, and dilution refrigerator. The cell is filled by condensing the 
desired mixture of helium isotopes through the fill tube. The 50 O resistor serves as the heater. 
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Fig. 2 - Numerical results for the relative change in the integrated ^He density integrated along 
the the neutron beam path as a function of input power P, parameterized as a/D; the lower plot 
shows the numerical derivative of the upper plot. These results are specific to the geometry of our 
experiment. 



Fig. 3 - Experimentally-determined relative change in scintillation rate data plotted as a function of 
X = aT^ . The data points for R{x)/R{0) > 0.6 represent a total count of about 30,000 scintillation 
pulses, from which a background of 1000 has been subtracted; the background is independent of ^He 
density and heater power. The errors in temperature and relative counting statistics are of order the 
size of the symbols used to plot the data. Solid lines indicate the slope at low heater power, while 
the dashed lines indicate the trends in the data at high power. 
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Fig. 4 - An expanded view of the low temperature data from Fig. |3[, which fall on a single line the 
the low power regime (circles, 0.45 K; diamonds, 0.55 K; squares, 0.65 K). This behavior indicates 
that D scales as T^^. A least-squares fit to a line y — mx + b gives b — 0.992 ± .004, m = 
(-4.77 ± 0.45) X 10~^s/cm^K^ with a reduced = 0.78. 



Fig. 5 - A comparison of our results for the mass diffusion coefficient D of ^He in superfluid ''He 
below 1 K with those of previous workers at higher temperature. 



